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Abstract 


The  collinear  crack  problem  in  an  inhomogeneous  orthotropic  medium  is 
considered  under  Mode  I  plane  strain  or  plane  stress  loading  conditions.  It  is  shown 
that  by  introducing  certain  averaged  orthotropy  parameters,  aside  from  a  scaling 
parameter  the  results  become  only  weakly  dependent  on  the  orthotropy  constants.  The 
main  results  of  the  study  consist  of  the  stress  intensity  factors  at  various  crack  tips  as 
influenced  by  the  material  inhomogeneity  parameter  and  by  the  relative  size  and 
position  of  the  secondary  cracks  with  respect  to  the  dominant  crack.  Also  considered 
is  the  crack/contact  problem  for  the  graded  medium  subjected  to  remote  tension  and 
bending  through  fixed  grips.  It  is  shown  that  the  crack  surface  contact  on  the 
compression  side  of  the  loading  has  a  magnifying  effect  on  the  stress  intensity  factor 
on  the  tension  side. 


1.  Introduction 

In  recent  years  there  has  been  considerable  interest  in  grading  the  thermomechanical 
properties  of  composites  as  a  new  tool  in  designing  materials  for  specific  applications  (  see 
Yamanouchi  et  al.  1990,  Holt  et  al.  1993  and  Ilschner  and  Cherradi  1995  for  review, 
applications  and  extensive  references).  Up  to  now  most  of  the  work  in  the  field  has  been 
on  metal/ceramic  particulate  composites.  By  proper  selection  of  the  constituent  materials, 
compositional  grading  and  the  processing  technique,  the  concept  may  be  used  to  develop 
new  materials  having  such  highly  desirable  and  seemingly  irreconcilable  properties  as  high 
heat,  corrosion  and  wear  resistance,  high  strength  and  high  toughness  in  the  same  material 
system.  The  technique  may  also  be  used  to  process  interfacial  regions  to  improve  bonding 
(Kurihara  et  al.  1990)  and  to  reduce  the  magnitude  of  residual  and  thermal  stresses  in 
bonded  dissimilar  materials  (Lee  and  Erdogan  1995,  Lee  and  Erdogan  1997). 

In  the  past  the  fracture  mechanics  studies  of  graded  materials  were  concerned 
primarily  with  the  influence  of  material  inhomogeneity  constants,  and  certain 
dimensionless  length  parameters  on  the  fracture  behavior  of  components  containing  a 
single  dominant  crack.  For  example,  the  model  and  the  mixed  mode  crack  problems  for 
an  infinite  isotropic  inhomogeneous  medium  were  studied  by  Delale  and  Erdogan  (1983) 
and  Konda  and  Erdogan  (1994),  respectively.  The  influence  of  the  length  parameters  on 
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the  stress  intensity  factors  and  the  strain  energy  release  rates  in  isotropic  graded  layers 
undergoing  spallation  and  surface  cracking  was  investigated  by  Chen  and  Erdogan  (1994) 
and  Erdogan  and  Wu  (1996),  respectively.  However,  because  of  the  nature  of  the 
techniques  used  in  processing,  the  graded  materials  are  seldom  isotropic.  Of  the  two  most 
commonly  used  processing  techniques,  generally  thermal  spray  would  give  a  lamellar 
(Sampath  et  al.  1995)  and  the  electron  beam  physical  vapor  deposition  a  columnar 
structure  (Kaysser  and  Ilschner  1995).  An  appropriate  model  for  such  graded  materials 
would  be  an  orthotropic  inhomogeneous  continuum  (Ozturk  and  Erdogan  1997).  One 
additional  factor  that  needs  to  be  considered  in  studying  the  basic  crack  problem  in  graded 
materials  is  the  so-called  multiple-site  cracking  or  the  crack  interaction.  In  homogeneous 
materials,  as  intuitively  expected,  the  crack  interaction  would  invariably  increase  the 
relative  magnitude  of  stress  intensity  factors.  In  graded  materials,  however,  such  a 
categorical  statement  is,  generally  not  possible  and  the  result  may  depend  on  the  relative 
positions  of  the  multiple  cracks  with  respect  to  the  direction  of  the  material  property 
variation  as  well  as  on  the  relative  dimensions. 

The  main  objective  of  this  article  is  to  study  the  influence  of  the  relative  size  and 
location  of  a  secondary  crack  on  the  stress  intensity  factors  in  an  orthotropic 
inhomogeneous  medium  containing  a  dominant  crack.  It  is  assumed  that  the  cracks  are 
collinear  and  are  located  in  a  plane  parallel  to  the  direction  of  material  property  variation. 
It  is  also  assumed  that  the  problem  is  one  of  plane  strain  or  generalized  plane  stress  and 
the  solution  of  the  elasticity  problem  in  the  absence  of  crack  is  known.  Thus,  the  plane  of 
the  cracks  is  one  of  symmetry  and  the  problem  is  a  model  crack  problem  in  which  the 
known  crack  surface  tractions  are  the  only  external  loads.  The  problem  is  formulated  for 
arbitrary  crack  surface  tractions  and  the  results  are  given  for  some  simple  loading 
conditions.  The  stress  intensity  factors  for  more  complex  loadings  may  then  be  obtained  by 
superposition.  Clearly,  the  result  would  be  valid  only  if  the  mode  I  stress  intensity  factors 
k\  at  all  crack  tips  are  positive.  A  negative  k\  implies  crack  closure  in  which  case  the 
problem  becomes  nonlinear  with  the  size  of  the  contact  zone  being  an  additional  unknown. 
This  problem  is  solved  by  assuming  that  the  graded  medium  contains  a  single  crack  and  is 
subjected  to  tension  and  bending  away  from  the  crack  region  through  rigid  grips  (  see 
Appendix  A  for  the  solution  in  a  homogeneous  isotropic  medium  which  is  obtained  in 
closed  form). 

2.  Formulation  of  the  Elasticity  Problem 

Let  Eii,  E22,  G 12  and  */y,  ( i,j  =  1, 2, 3)  be  the  engineering  elastic  parameters  for  an 
orthotropic  inhomogeneous  plane  x\ ,  X2.  To  replace  them  we  introduce  the  following  four 
parameters  (Krenk  1979,  Cinar  and  Erdogan  1983): 

E  =  VES,  <'=VwTu  =  I11  =  — ,  k=-£--v  (1) 

&22  ^21  ^Yl 

for  plane  stress  and 
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_ •£'11-^22 _ 

(1  —  1^13^31  )(1  —  ^23^32) 


(l/j2  +  1^13  ^32) (^21  +  ^23  ^31 ) 

(1  —  V\z  ^3l)(l  —  ^23  ^32 )  /  ’ 


,4  _  Exx  (1  ~  ^23^32)  _  E  _ 

°  E22  (1  -  V13V3iy  K  2 g12  *  C) 

for  plane  strain  conditions.  In  graded  materials  generally  E,  v,  k  and  5o  would  be 
functions  of  xx,x2.  However,  previously  it  was  shown  that  the  results  are  only  weakly 
dependent  on  v  (Delate  and  Erdogan  1983,  Ozturk  and  Erdogan  1997)  and  it  may  be 
assumed  that  Exx,  E2 2  and  GX2  vary  proportionately.  Consequently,  the  parameters  k  and 
60  may  be  assumed  to  be  constant  and  the  material  inhomogeneity  may  be  represented  by 
the  function  E(x x,x2)  only.  Note  that  60  =  1  and  k  =  1  correspond  to  the  isotropic 
medium  and  for  k  <  —  1  the  problem  has  no  feasible  solution  (Erdogan  and  Wu  1993).  By 
introducing  6q  as  a  scaling  parameter,  in  the  usual  notation  we  define 

x  =  y  =  yfc  x2, 

u(x,y)  =  u\  {xi,x2)y/6o,  v(x,y)  =  U2{xX,X2)/ y/~6 0, 

C Txx(x,y )  =  an(xi,x2)/60,  ayy{x,y)  =  60a22(xx,x2), 

<Txy(xty)  =  <T12(xx,x2).  (3) 

Furthermore,  by  assuming  that  the  stiffness  E  varies  in  xi  direction  only  and  in  the  crack 
region  may  be  approximated  by  (Fig.  1) 

E{x  x,x2)  =  E{x  1)  =  EoeQXl  =  E(x)  =  Ebe7*,  7  =  ay/fo,  (4) 

the  equilibrium  equations  may  be  expressed  as 


d2u  d2u  d2v 
dy 2  1  dx2  dxdy 


d2v  0  d2v  0  d2u 

d*+PlW  1 


( du  dv\ 

Pn{di+,,di)=0’ 

(du  dv\ 

1{d^  +  di)  ”0’ 


where 


2  (k  +  v) 


2=1+  V@1  • 


Due  to  symmetry  by  considering  y  >  0  half  plane  only,  from  (5)  it  may  be  shown  that 


1  r°°  4 

«(*, : V)  =  W  „  E  eX)*~ikx  dk, 

1  r°°  4 

v{x,  y)  =  eXjV~ikx  dk,  0  <  y  <  00, 
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«•=(*.*)  =  -  iAjk)Bj  ex’y~ikx  dk, 

*„(*. »)  =  ^  ^/“  I>>  -  ^~ikx  dk • 

<r„(*.  v)  =  2(f +t)  ~  ik)Bi  eX‘’~itZ  dk’0<V<°°-  <8> 

where 

Ai  =  -  A3  =  |^(C2  +  2*77)  +  ^\/(C2  +  2/«r?)2  -  4t;2|  , 

A2  =  -  A4  =  1 1  (C2  +  2kt;)  -  i^/(C2  +  2kt?)2  -4^j  ,  (9) 

C2  =  I'T2,  77  =  A;2  +  nfc,  (10) 

and 


a  (*&*  “  7/5i^)  , 

j  ~  A2  -  far,  Xj 


3  =  3,4. 


(11) 


The  unknown  functions  S3  and  S4  are  determined  from  the  following  boundary 
conditions 


^12(^1, 0)  =  0,  —  00  <  X\  <  00.  (12) 

n 

Ototei.+o)  =  pM,  Xi  6  L,  L  =  Y,Lk,  Lk  -  (akbk), 

1 

W2(^i>+0)  =  >/i5o v(x, +0)  =  0,  xi  e  L',  L  +  L'  =  (- 00,00).  (13) 

Here,  it  is  assumed  that  the  crack  surface  traction  p(x  1)  is  a  known  function  and  the 
medium  contains  n  collinear  cracks  along  x2  =  0  (Fig.  1). 

By  defining  the  following  new  unknown  function 

Q 

— u2(xi,+0)  =y(xi),  -oo<zi<oo,  (14) 

UX 1 

it  can  be  shown  that  the  boundary  conditions  reduce  to  (Ozturk  and  Erdogan  1997) 

-  /  [7—^ —  +  S2(ti  -  xi)  g(ti)dti  =  •■Sl  e~axi p{x  1 ),  xi  €  L,  (15) 

7r  JL  Lh  ~  xi  J  So 
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where 


-/  g{xi)dxi=0,  j  =  l,...,n  (16) 

njLj 

Si  =  y/n  +  Kl,  S2  =  y/n  -  «1,  Kl  =  V«2  -  1,  (17) 

Ri{t\  —  %i)  ~  2  /^~^1  — 

i?l(s-x)  =  f  2  [3ft(h(A;))cos(fc(s  —  x))  —  9‘(/i(fc))sin(A:(s  —  x))j  dfc,  (18) 


.[(sH-s 

*L  Ax 


s2)(fc  +  «7) 

i  +  a2 


/i(fc)  = 


2(fc  +  fry) 

V72  +  4(/c2  +  *7/:) 


/c  7^  1 


1  ,  «  =  1. 


In  real  materials  —  1  <  k  <  coand  si  +  s2  is  always  real.  Also,  for  7  =  0,  Ai  =  ks\, 
A2  =  /cs2,  h(k)  =  0,  R-2  =  0  and  (15)  would  reduce  to  the  integral  equation  for  a 
homogeneous  orthotropic  medium. 

It  is  interesting  to  note  that  in  the  special  case  of  isotropic  inhomogeneous  materials 
the  integral  in  (18)  can  be  evaluated  in  closed  form  giving  (Appendix  B) 

Mh  -*i)  =  ^|i'o!|/fi(i'ol?l)  +  /fo(^oKI)}  - 

«-§&-».).  (20) 

where  K0  and  K\  are  modified  Bessel  functions.  Also  note  that  the  kernel  i?2  depends  on 
the  inhomogeneity  constant  a  and  the  elastic  constants  v  and  k  but  not  on  E$  and  6q. 
Similarly  si  and  s2  depend  on  the  shear  parameter  k  only.  It  may  then  be  concluded  that  in 
(15)  Eo  and  60  are  simply  scaling  constants  and  would  have  no  influence  on  the  stress 
intensity  factors  which  may  be  defined  by  and,  after  solving  (15),  evaluated  from 

ki(bk)=  lim  v/2(xj  -  6*) ^22(^1, +0) 

xi  — >  ojt  +0 

-  -j 3{Xi)- 

h(ak)=  lim  V2(a*  -  Zi)<722(a:i,+0) 

xi  -*  ak  -0 

_  E(x\)  fKTZ  ~  \  l _ 1  ~  /on 


Iim  n7  T  U  =  (21) 

xi^>ak+0{Si  +  S2)Oo 
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(22) 


Also,  from  (13)  and  (14)  the  crack  surface  displacements  may  be  obtained  as 

rxi 

=  /  9$i)dt\,  dj  X\  <!  bj,  j  ==  1,...,ti. 

Jdj 


3.  On  The  Solution  Of  The  Integral  Equation 

Since  the  closed  form  solution  is  not  available  the  integral  equation  (15)  must  be 
solved  numerically.  A  simple  and  very  effective  way  to  do  this  would  be  to  reduce  (15)  to 
a  system  of  n  integral  equations  each  having  the  support  (—1,  1)  by  defining 

2ti  —  ( bk  +  Ofc) 

q  =  — -  —j — ,  ak  <  ti  <  bk,  -l<q<l, 


2xi  -  { bk  +  ak ) 


ak  <  xi  <  bk,  -1  <  r  <  1, 


(6* -a*)  ’  — 

<t>k{g)  =  9(h),  ak  <  ti  <  bk,  -1  <  q  <  1, 


fk(r)  =  p(xi)/P0, 


ak  <  x i  <  bk,  -1  <  r  <  1, 


RikMq  -  zik ))  =  i?2(<i  -  xx)/ck,  a,i  <x i  <  bi,  dk  <  h  <  bk, 


—  1  <  9  <  1» 


Ck  =  ~  a*)>  ^ k  =  ^(^k  +  a*), 


Zik  =  ~ r  +  —  (di  -  dk),  (i,k)  =  -1  <  r  <  1,  (23) 

Ck  Ck 

where  P0  which  has  the  dimension  of  stress  is  a  normalization  constant.  The  equations 
(15)  and  (16)  may  then  be  expressed  as 


i  r  f-L 

nJ- 1  L9- 


+  CiRii\<f>i(q)dq 


+  Y f  f — - \-ckRik  <t>k{q)dq 

is' nJ-i[q- Zik  J 


(si  +  S^SqPq 
E0exp(a(cir  +  dt-)) 


fi(r), 


-1  <  r  <  1,  i  =  1 


J^(f>k{q)dq  =  0,  k  — 


The  system  of  integral  equations  may  be  solved  in  a  standard  manner  by  expressing 
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<f>M)  = 


(si  +  52)^0-Pq 

E0exp(adk) 


[Z - O  'y  1  AkmTm{ci), 

Vl-?2^o 


— 1  <  q  <  1,  k  =  n, 


(26) 


where  Tm(q )  is  the  Chebyshev  polynomial  of  the  first  kind  and  Akm  are  unknown 
dimensionless  constants.  From  (25)  and  othogonality  conditions  of  Tm's  it  follows  that 
Ak o  =  0.  The  coefficients  Akm  are  determined  by  substituting  from  (26)  into  (24)  and  by 
using  the  method  of  reduction  (Kantorovich  and  Krylov  1958,  Erdogan  1978).  The  stress 
intensity  factors  may  then  be  obtained  from  (21)  and  (26)  as  follows: 

OO 

h(h)  =  A*-, 

m=l 

oo 

*i(o*)=  Poyft~ke-aCkY,(-l)mAkm,  k  =  (27a,  b) 

m= 1 


4.  Results  and  Discussion 

The  main  interest  in  this  study  is  in  evaluating  the  influence  of  relative  position,  size 
and  distance  of  a  "small"  crack  lying  in  the  plane  of  a  dominant  crack  on  the  stress 
intensity  factors  in  an  inhomogeneous  orthotropic  medium.  With  the  additional  knowledge 
of  subcritical  crack  growth  characteristics  of  the  material,  one  may  then  be  able  to 
determine,  for  example,  fatigue  or  corrosion  crack  growth  rate  at  each  crack  tip.  The 
problem  is  solved  under  two  sets  of  loading  conditions.  In  the  first,  it  is  assumed  that  the 
crack  surface  tractions  are  given  by  (see  13) 

p(®i) =  Po  Pi  (”■)  P2  (“■)  ■"»(“)»  XieL,  (28) 

where  Po,PuP2  and  p3  are  the  measure  of  the  magnitude  of  external  loads  and  are  known 
constants  and  a  is  a  normalizing  length  parameter  (usually  the  half  length  of  the  dominant 
crack  ).  Since  the  problem  is  linear,  these  results  may  be  helpful  to  find  an  approximate 
solution  to  the  crack  problem  for  a  given  specific  loading.  The  second  external  load 
considered  is  a  remote  displacement  loading  of  the  form 

€22  (xi ,  x2 )  =  e0  +  €i  (xi/a)  (29) 

giving 

P(x i)  =  -  ^-eaXl  (e0  +  €i  (~))>  xiEL.  (30) 

Summary  of  the  results  for  an  isotropic  inhomogeneous  medium  containing  a  single 
crack  along  x?,  =  0,  —  a  <  Xi  <  a  is  shown  in  Tables  1-3.  These  results  are  both  accurate 
and  comprehensive.  Tables  4-6  show  the  comparison  of  the  stress  intensity  factors 
obtained  in  this  study  with  that  found  by  Delale  and  Erdogan  (1983).  From  the  solution 
given  in  the  previous  section  the  crack  surface  displacement  may  easily  be  expressed  as 
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w2(zi,+0) 


(31) 


where  u2(xi,+0)  =  v(x,0), 

v0  =  2aPo(l  —  v2)/Eq  (32) 

for  isotropic  and 

Vo  =  (5i  +  s2)60P0a/ Eq  (33) 

for  orthotropic  inhomogeneous  materials  (Eq  =  Eq  for  plane  stress  and 
Eo  —  Eo/(l  —  w2)  for  plane  strain  ).  Some  sample  results  for  the  crack  surface 
displacement  are  shown  in  Figure  2. 

The  calculated  plane  strain  results  for  the  normalized  stress  intensity  factors  at  various 
crack  tips  in  an  isotropic  inhomogeneous  medium  containing  two  unequal  cracks  under 
fixed  grip  (e22  =  e0)  or  fixed  load  (cr^xi,  0)  =  — po)  conditions  are  given  in  Figures  3- 
16.  Some  of  the  limiting  cases  of  these  results  for  (ci  +  c2)/d  — >  0  (the  uncoupled  case) 
and  (ci  +  c2)/d  — >  1  (  the  case  of  a  single  crack  of  length  2(ci  +  c2))  may  be  found  in 
Table  7.  One  may  note  that  as  (ci  +  c2)/d  — *  1  or  a2  — ►  &i,  k\(bi)  — >  oo,  k\(a 2)  —*  00 
and  &i(ai)  and  ki  (b2 )  approach  the  single  crack  values.  As  shown  by  Boduroglu  and 
Erdogan  (1983)  and  as  may  be  observed  in  Figure  7,  for  02  —*  b\  fci(oi)  and  k\  (b2 )  are 
highly  ill-defined  ( that  is,  their  slopes  are  unbounded  as  (ci  +  c2)/d  — ►  1). 

Some  sample  results  for  stress  intensity  factors  in  an  orthotropic  inhomogeneous 
medium  containing  two  unequal  collinear  cracks  under  plane  strain  conditions  are  given  in 
Figures  17-21  and  Tables  8  and  9.  Aside  from  the  geometric  parameters  c2/ci  and 
(ci  +  c2)/d,  as  pointed  out  previously,  the  main  variables  in  this  case  are  the  material 
inhomogeneity  constant  a  and  the  shear  parameter  k.  Both  the  figures  and  the  tables  show 
that  the  dependence  of  the  stress  intensity  factors  on  a  is  strong  and  on  k  is  rather  weak. 

The  results  for  the  crack/contact  problem  under  remote  loading  e22  —  £0  +  £1  ix/a) 
are  shown  in  Tables  10,  11  and  Figures  22-25.  Figure  22  shows  the  crack  surface 
displacement  for  £0/ei  =(3  =  0  (remote  bending)  obtained  by  ignoring  the  crack  surface 
interpenetration.  Thus,  on  the  compressive  side  of  the  external  load  the  solution  gives 
negative  displacement  which  is  physically  unacceptable.  Figures  23  and  24  show  the  crack 
surface  displacement  for  (3  =  0  obtained  by  taking  into  account  the  smooth  contact  of  the 
surfaces.  Note  the  slight  dependence  of  the  crack  surface  displacement  on  the  sign  of  e\. 
Table  10  shows  the  stress  intensity  factors  and  the  size  of  the  contact  zone  again  for 
remote  loading  e22  =£i(xi/a)  and  for  various  values  of  the  material  inhomogeneity 
constant  a,  (E(xi)  =  E0exp(axi),  see  (4)).  The  contact  occurs  along  x2  =  0, 
-a  <  x  1  <  —  b,  (-a  <  -b)  for  ^  >  0  and  along  x2  =  0,  b  <  x\  <  a,  (b  <  a)  for 
£1  <  0.  The  table  also  shows  the  (theoretical)  stress  intensity  factors  calculated  by 
disregarding  the  crack  surface  penetration  (no  contact  case).  Note  that  the  presence  of 
contact  tends  to  magnify  the  stress  intensity  factor  at  the  crack  tip  on  the  tension  side  ( 
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k\(a)  for  £\  >  0  and  ki(—a)  for  £1  <  0).  Note  also  the  highly  significant  influence  of  the 
material  inhomogeneity  constant  on  the  stress  intensity  factors. 

Some  sample  results  for  the  combined  loading  "tension"  and  "bending", 
s22  =  e0  +  e1(xi/a),  (eo/ei  =  /?)  are  shown  in  Figure  25  and  Table  11.  Figure  25  shows 
the  crack  surface  displacement  for  various  values  of  j3.  For  (3  >  (3C  and  the  contact  zone 
size  becomes  zero,  that  is  b  —  a,  ft  =  (3C  corresponding  to  k\{—a)  =  0.  Also,  b  <  a  for 
(3  <  Pc- Table  11  shows  the  crack  contact  zone  sizes  and  the  corresponding  stress 
intensity  factors.  If  the  material  is  homogeneous  and  isotropic  the  crack/contact  problem 
can  be  solved  in  closed  form  which  is  given  in  Appendix  A. 

Finally,  a  sample  result  giving  the  distribution  of  ayy(x,0)  along  the  net  ligament 
b\  <  x  <  02  is  shown  in  Fig.  26.  Note  that  the  solution  given  in  this  study  is  that  of  the 
perturbation  problem  which  must  be  added  to  the  results  obtained  from  the  uncracked 
medium  in  order  to  find  the  total  solution. 

5.  Conclusions 

In  the  case  of  two  collinear  cracks  with  crack  tips  at  (01,61)  and  (02,62)  and 
ai  <  61  <  a2  <  62,  61  -  ai  »  62  -  02  (that  is,  (01, 61)  being  the  dominant  crack ),  if  the 
medium  is  homogeneous  (isotropic  or  orthotropic)  then  under  mode  I  loading  generally 
k\{b\)  >  k\ (oL)  >  £1(02)  >  ki(b2),  meaning  that  fastest  crack  growth  would  take  place 
at  the  inner  tip  of  the  dominant  crack.  As  the  ligament  a2  -  61  decreases  and  tends  to  zero 
ki(a2)  as  well  as  ki  (61)  becomes  unbounded.  However,  if  the  medium  is  inhomogeneous, 
the  values  of  the  stress  intensity  factors  may  be  influenced  quite  significantly  by  the 
relative  positions  of  the  dominant  and  the  secondary  cracks  with  respect  to  the  direction  of 
material  property  variation. 

After  introducing  the  "averaged"  orthotropy  parameters  E,  u,  k,  and  So,  in  mode  I 
crack  problems  for  an  inhomogeneous  orthotropic  medium  E(0)  and  <50  act  as  scaling 
constants  and  v  and  k  do  not  seem  to  have  a  significant  influence  on  the  stress  intensity 
factors.  In  crack  problems  involving  graded  materials  the  dominant  factor  remains  to  be 
the  material  inhomogeneity  constant. 

In  mode  I  and  mode  HI  crack  problems  for  isotropic  graded  materials  described  by 
fi(x)  =  /x0exp(ax)  the  kernels  of  the  integral  equations  can  be  evaluated  in  closed  form. 
Consequently,  for  these  materials  highly  accurate  benchmark  solutions  can  be  obtained.  In 
this  article  the  mode  I  stress  intensity  factors  in  such  materials  are  provided  for  tension  and 
bending  through  fixed  grips  as  well  as  third  degree  polynomial  tractions  acting  on  the 
crack  surfaces. 
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APPENDIX  A 


Crack  Closure  for  a  Homogeneous  Isotropic  Medium 


Let  the  medium  contain  a  single  crack  of  length  2 a  along  the  x  axis  and  be  subjected 
to  crack  surface  tractions  axy(x,  0)  =  0,  ayy(x,0)  =  p{x),  —a<x<a.  By  defining 
dv(x,  +0)/dx  =  g(x),  the  integral  equation  of  the  problem  may  be  expressed  as 


1  fa  g(s)ds  2 

Z  /  "7 — r*  TtP(x)>  ~a<x  <a, 

7T  J —a  S  X  £/q 


subject  to 


^  f  g{s)ds  =  0, 
7T  J —a 


where  Eq  =  E  for  plane  stress  and  Eq  =  E/(  1  —  u2)  for  plane  strain  conditions.  By 
substituting 

x  =  ax' ,  s  —  as',  g(s )  =  G(s'),  p(x)  =  P(x').  A3 

A1  and  A2  become 

l  f  =  J-P(A  f  G(s’)  ds‘  =  0.  A4 

7T  j—l  ^  X  -tl/O  •/— 1 

We  now  assume  that  the  solution  of  (A4)  is  of  the  form 


G(s’)  =  £ 


AnTn(s') 


o  V  l  -  s' 


Considering  the  orthogonality  condition 


1,  m  =  n  =  0 

1  f  Tn(q)Tm(q)  ,  ,1/0  _  \  -I 

~  — /,  0  dq  =  <  1/2,  m  =  n>i, 

ft  J- 1  >/l  - 

0,  m  ^  n, 


from  (A4b)  it  may  be  seen  that  Ao  =  0.  Also,  by  substituting  from  (A5)  into  (A4a)  and  by 
using  the  properties 


1  f1  _ Tn 

nJ-i  ( q-r )- 


(q-rj^/l-g2 


0,  n  =  0,  |r|  <  1 

dq  =  <  Un-i{r),  n  >  1,  |r|  <  1, 

Gn{r),  n  =  0, 1, ... ,  \r\  >  1, 


12 


Tn(q)  =  cos (nd),  Un(q) 


sin  (n  +  1)0 


cos0  =  q, 


Gn  (r)  = 


(r  ~  (|r|/r)v^T) 


we  obtain 


YjAmUm-l{x')  =  2rP{x').  A1 

m=l  -®0 

As  an  example  consider  now  the  remote  "bending"  of  the  medium  for  which  we  have 
P{x,)  =  -E0e1x' = -E0eiUi(x')/2.  A1 

From  (A10)  it  then  follows  that 

A\  =0,  A2  =  —£1,  Am  =  0,  m  >  3.  A1 

and  the  solution  becomes 


G(s') 


£1^2  (s')  £i(2s'2  —  1) 


1  -  s' 


giving  the  stress  intensity  factors 


h(a)  =  ^E0£iy/a,  ki(-a)  =  ~^E0£iy/a.  A14 

Since  the  negative  stress  intensity  factor  (A14b)  implies  interpenetration  of  the  crack 
surfaces,  the  solution  given  by  (A13)  is  not  valid.  The  correct  solution  may  be  obtained  by 
taking  into  account  the  crack  surface  contact  near  the  end  x  =  —a.  Let  the  contact  region 
be  —a<x<—b  with  b  being  an  unknown  constant  to  be  determined  from  the 
smoothness  condition  ki(—b)  =  0.  The  integral  equation  (Al)  and  the  single-valuedness 
condition  (A2)  may  then  be  modified  as 


l  r  g(8)d8  =  2^  -b<x<a, 

7T  J —})  S  X  Eq 

1  r 

-  g(s)  ds  =  0, 

7t  J-b 


By  defining  the  following  normalizations 

.  _  _  _  ct  +  6, 

s  =  cS  +  d,x  =  cX  +  a,  c  =  — - — ,d  - 

u 


a  —  b 


,g(s)  =  G(S),p(x)  =  P(X),  A17 


(A15)  becomes 
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i£mds=2_P{X)'  _1<x<1 


S-X  E0 ' 


Applying  again  remote  "bending"  we  have 


p{x)  =  -E0£i(-)  =  -E0e  1  [-X+  -1  =  P(X) 
\a/  [a  aj 


The  integral  equation  (A18)  then  becomes 
1  flG(S)dS  fc  d] 

«J-l-S=T  =  -2£'[ZX+a\’ 


•1  <  X  <  1, 


Again  assuming  the  solution  as  given  by  (A5)  it  can  be  shown  that 


OO 

Y)  =  -2e,  -JC+-  . 

L»  “J 


By  observing  that  t/0(I)  =  1  and  Ux  (X)  =  2X ,  from  (A21)  it  follows  that 

Ai  =  —  2ei-,  A2  =  —  £1-,  Am  —  0,  m  >  2.  A22 

a  o 

Thus,  the  solution  becomes 

G(s)“-a7r=^[MS+c(1“2SS)]' 

giving 

h  (a)  =  E0ei  x/c  ^  M ,  A24 

*1  (-6)  =  -E0ei  y/c  ^  •  A25 

From  the  smooth  contact  condition  ki  (—6)  =  0,  the  unknown  constant  b  is  found  to  be 
,  1 

0  =  go  A26 

From  A24  and  A26  it  then  follows  that 

k\  (u)  2  k\  ( cl )  2  /  . 

PoeiV^  =  3  °r  P0£i\A  =  3^^  =  °-5443,  A27 

Figure  A1  shows  the  crack  surface  displacements  obtained  from  (A13)  by  ignoring  the 
interpenetration  of  crack  surfaces  and  from  (A23)  by  taking  into  account  crack  closure 


along  -a  <  x  <  -b.  From  Figure  A1  and  equation  (A27)  it  is  seen  that  the  crack  closure 
near  the  end  x  —  -a  leads  to  magnification  of  stress  state  near  the  crack  tip  x  =  a  and 
there  is  approximately  a  9%  increase  in  the  stress  intensity  factor  k\  (a)  over  the  nominal 
value  obtained  by  ignoring  crack  surface  contact. 


x/a 


Figure  A1  Crack  surface  displacements  in  a  homogeneous  isotropic 
medium  subjected  to  remote  bending  eyy  =  £i(x/a)  and  obtained  by 
ignoring  crack  closure  and  by  taking  it  into  account;  vq  =  2e\a. 


After  solving  the  perturbation  problem  (A15)  the  stress  (Tyy(x,  0)  on  the  plane  of  the 
crack  may  be  obtained  by  using  (A7),  (A9),  (All)  and  (A18)  and  by  adding  the 
homogeneous  stress  ayy  =  E0Si(x/a)  as  follows: 


(*yy{xi  0) 


0,  -b  <  x  <  a, 


^[A1Gl(X)  +  A2G-l(X)}  +  Eae1 


X  = 


x  -  d 


c  = 


a  +  b 
2  ’ 


—  oo  <  x  <  —b,  a  <  x  <  oo,  A28 


or  by  observing  that  b  =  a/3  it  may  be  shown  that 
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<%(*>  0) 
Eq£\ 


2X  + 1  IV1  , 

+  —3—.  \X\  >  1.  A29 

Figure  A2  shows  the  stress  distribution  expressed  by  (A29).  Note  that,  as  expected, 
a-yy(x ,  0)  <  0  for  x  <  -b. 


Figure  A2  The  stress  ayy  in  the  plane  of  the  crack  in  a  homogeneous 
isotropic  medium  under  remote  "bending"  eyy  =  £1  {x/a). 


Let  us  now  assume  that  a  homogeneous  isotropic  medium  containing  a  crack  along 
-a  <  x  <  a,  y  =  0  is  subjected  to  remote  “tension”  eyy  —  e0  as  well  as  “bending” 
£yy  —  £i  (x/a).  The  input  function  in  the  integral  equation  (Al)  would  then  be 


p{x)  —  —Eo  jfo  +  £1  (“)]  —  ~Eq£i  [/?  + 

(:)]• 

A30 

giving  the  solution  as  follows: 

G(s')=  A-^ArTiM+AiTiMl 

s'  =  s/a, 

A31 

Ai  =  -2e1f3,  A2  =  -Ei, 

A32 
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A33 


Figure  A3  shows  the  crack  surface  displacement  v(x,  0),  —b<x<a,  the  size  of  the 
contact  zone  b  and  the  normalized  stress  intensity  factors  for  various  values  of  the  strain 
ratio  /?  =  eo/ei. 
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Figure  A3  The  crack  surface  displacement,  the  contact  zone  size  and  the 
normalized  stress  intensity  factors  in  a  cracked  isotropic  homogeneous 
medium  subjected  to  the  remote  loading  eyy  —  e0  +  £\{x/a),  v0  =  2e\a. 
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APPENDIX  B 


Evaluation  of  the  Kernel 


For  an  isotropic  inhomogeneous  medium  containing  a  series  of  collinear  cracks  along  the 
x  axis  the  kernel  of  the  integral  equation  (15)  may  be  expressed  as 


-  !£ 


k  +  ia 

'voi1  +  A(k2  +  iak) 


eik“dk 


=  2(- 


+  i?2M),  U)  =  s  -  X. 


Defining  A  =  a/2  and  observing  that  k2  +  iak  —  (k  +  iA)2  +  A2,  from  (Bl)  it  follows 
that 

D,,  JT  (k  +  i2A) 


r^) =~  r  -  (fc+t2A)  . — t^dk 

lJ~ 00  J (y  +  \)A2  +  (k  +  iA)2 


Changing  the  variables  k  =  AV  and  assuming  A  >  0,  (B2)  becomes 

R(u)  =  -  [°°  -  +  H  eiAV“dV, 

1  J-oc  J(u  +  l)  +  (y  +  i) 2 


or,  by  letting  U  =  V  +  i,  u  +  1  =C2  and  B  —  Au>  we  find 


00+4  u  +  i 

oo+i  y/CPTU2 


eBUdu. 


By  substituting  U  =  —U'  and  observing  that 


izEL  ±  il  Q~iBU‘ 

/C2  +  U'2 


/OC 

. 


*  ~  i)  Q-iBU' 

\/C2  +  U'2 


(B4)  may  be  written  as 


R(u> )  =  Ae£ 


=  2  Aee9fc 


i  V  *v 

re 

(  /*oo+t 

1/  i 


(g  - 

>C2  +  u2 

(£±j) 
Jc‘ 2  +  u2 


eBUdu 


gjH) 

iJ&VU2 


elBUdU  >. 
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We  now  consider 


roo+i  (U  +  i) 


iJC2  +  U2 


tiBUdU,  B>  0. 


;(R+1) 


i(C  +  e) 
iC 
KC  -  e) 


R+  i 


Figure  B1  Contour  used  in  evaluating  Ii 


Referring  to  Figure  B1  we  have 


(z  +  i) 
/z2+C2 


p R+i  p  pi(C+e)  p  pi 

tiBzdz  =  +  +  +/  +  /  =0> 

Ji  Jr  Ji(R+ 1)  J  7  J  i(C—e) 


For  R-*  oo  and  e-»0  it  can  be  shown  that  the  integrals  over  the  arcs  7  and  T  in  (B8)  are 
zero  and  from  (B8)  it  follows  that 


iJC2  +  U2 


Figure  B2  Contour  used  in  evaluating  I2 
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Let  us  now  assume  that  B  <  0,  B  =  —B'  and  define  the  function 


roo+i  (U  +  i) 


Ji  Xy/C 2  +  IP 

Referring  to  Figure  B2  we  have 


(*  +  *) 
/C2+Z2 


e~lBzdz 


Q~iBfudU,  B'  >  0, 


p± t+t  p  p— i(C+e)  p  pi 

=  +  /  +  /  +  /  +  /  =°- 
7r  J-HR-1)  A  J-i(C-e) 


(R- 1)  ^7  ./  —i(C—e) 


Again,  for  i?-+oo  and  e-^0  it  can  be  shown  that  the  integrals  over  the  arcs  7  and  T  in 
B 15  are  zero  and  I2  becomes 


f  +  ’2  e-igudU 

iy/C' 2  +  U2 


-I'. 

-j. 


(*±0  -i&z 


iC  iy/C2  +  z2 
*  (*  +  *) 
-iC  iyf C2  +  z1 


e  Zdz 


e~tBzdz 


~L  v£=we~*’dv  8,6 


(U  +  i)  .„«  I 

iyfc 2  +  C/2  f 


72  -C2 


1  =  r  (.~y  +  1}  e-^dy 
J  7c  Vy2  -  <?2 

_  r°  (-cY  +  i)  gcy 

~A 


=  -CKxiB'C)  + 


Thus,  for  any  real  constant  positive  or  negative,  we  find 
J*(w)  =  2Acb^C^K1(\B\C)  +  ^0(|B|C')|  =  2^  +  R2(uj)\ 

A  =  a/2,  C  =  y/  l  +  i/,  £  =  Au;  =  ^(s  —  :r). 

A 
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APPENDIX  C 


The  Antiplane  Shear  Problem  for  a  Graded  Medium  Containing  a 

Crack 


It  can  be  shown  that  the  integral  equation  corresponding  to  an  infinite  inhomogeneous 
isotropic  medium  containing  a  crack  and  subjected  to  anti-plane  shear  loading  is  given  by 


-  f  — —  +  Ma(x,  s)  g(s)  ds  =  — e  ax<ryz(x,  0), 
7T  J-a  Ls  ~  x  J  A*0 


—a<x<a, 


Cl 


9(x)  = 


dw(x,  0) 

di  ’ 


Ma(x,s)  =  fee|J|i^i(|B|)  +  K0(\B\)}  ~ 


«=§(*-*) 


C2 

C3 

C4 


where  Kq{z)  and  K\{z)  are  the  modified  Bessel  functions  of  order  zero  and  one 
respectively.  It  is  again  assumed  that  the  shear  modulus  of  the  medium  is  approximated  by 
fi(x)  =  /i0exp  (ax).  Note  that  the  kernel  of  the  integral  equation  (C3)  is  nearly  identical  to 
that  given  by  (20)  for  the  in-plane  problem.  Table  Cl  shows  some  sample  results  obtained 
from  (Cl)  and  their  comparison  with  that  obtained  by  delale  (1985)  numerically.  Since  the 
kernel  is  in  closed  form,  (Cl)  can  be  solved  within  any  desired  degree  of  accuracy. 


Table  Cl  The  normalized  mode  III  stress  intensity  factor  for  an 
inhomogeneous  isotropic  medium  subjected  to  uniform  crack  surface 
traction,  ay2(x,  0)  =  —  po,  -a  <  x  <  a. 


Delale  [1985] 

Present 

aa 

k3(a) 

fc3(-o) 

k3(a) 

h(-a) 

Poy/a 

Poy/a 

Poy/a 

Poy/a 

0.0 

1.0 

1.0 

1.0 

1.0 

0.1 

1.024 

0.973 

1.0228 

0.9731 

0.2 

1.045 

0.944 

1.0427 

0.9443 

0.3 

1.063 

0.914 

1.0605 

0.9149 

0.4 

1.080 

0.884 

1.0763 

0.8857 

0.5 

1.095 

0.855 

1.0906 

0.8570 

0.75 

1.127 

0.787 

1.1205 

0.7894 

1.00 

1.153 

0.726 

1.1438 

0.7291 
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Figure  1  Geometry  and  notation  of  the  collinear  crack  problem 
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Figure  2  Crack  surface  displacements  v(x,0)  in  an  inhomogeneous  isotropic  medium 
under  uniform  pressure  po  applied  to  the  crack  surfaces.  (  v0  =  2apo(l  —  v2)/Eq, 
v  =  0.3,  plane  strain  conditions.) 
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5 


(c1+c2)/d 


Figure  3  Stress  intensity  factors  at  the  tips  ai  and  fciof  two  unequal  cracks  in  an  isotropic 
FGM  under  fixed-grip  condition,  ayy(x,0)  = —Eoeoexp(ax),  ko  =  Eoeoy/ci, 
E0  =  E0/{1  -  v2),  c2/Cl  =  0.25,  i/  =  0.3. 


Figure  4  Stress  intensity  factors  at  the  tips  oi  and  &i  of  two  unequal  cracks  in  an  isotropic 
FGM  under  fixed-grip  condition,  ayy(x,0)  = -E0eoexp(ax),  k0  =  Eoeoy/ci, 
E0  =  Eo/(l  ~  v\  (c2  +  cx)/d  =  0.75,  i/  =  0.3. 
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(ci  +c2)/d 

Figure  5  Stress  intensity  factors  at  the  tips  ai  and  b\  of  two  unequal  cracks  in  an  isotropic 
FGM  under  fixed-grip  condition,  ayy(x,0)  =  -£70e0exp(ax),  ko  =  Eoeoy/ci, 
E0  =  E0/(l  -  u2),  c2/ci  =  0.25,  v  =  0.3. 


(ci  -\-c2)/d 


Figure  6  Stress  intensity  factors  at  the  tips  ai  and  b\  of  two  unequal  cracks  in  an  isotropic 
FGM  under  fixed-grip  condition,  ayy(x,0)  = -Eo£oexp(ax),  ko  =  Eo£oy/ci, 
Eq  =  Eo/(  1  —  v2),  aci  =  0.5,  v  =  0.3. 


26 


Figure  7  Stress  intensity  factors  at  the  tip  a\  of  two  unequal  cracks  in  an  isotropic  FGM 
under  fixed-grip  condition,  ayy(x,0)  = -Eoeoexp(ax),  ko  =  Eoeo^/ci, 
Eo  =  Eo/(l  -  v2),  aci  =  0.5,  v  =  0.3. 


ac\ 

Figure  8  Stress  intensity  factors  at  the  tips  02  and  b2  of  two  unequal  cracks  in  an  isotropic 
FGM  under  fixed-grip  condition,  <ryy(x,0)  =  -Eo£oQ\p(ax),  ko  =  Eoeoy/ci, 
E0  =  JEb/(l  -  v2),  c2/c  1  =  0.1, 1/  =  0.3. 
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Figure  9  Stress  intensity  factors  at  the  tips  ai  and  &iof  two  unequal  cracks  in  an  isotropic 
FGM  under  fixed-grip  condition,  cryy(x,0)  =  —Eo£o^p(o:x),  k0  —  Eq£q  yfc[, 
Eo  =  E0/(l  —  v2),  cz/ci  =  0.25,  v  =  0.3. 


MM/*  o 


c2/c,=  1 .0 
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0.25 

0.0 


M°i)/fco 


Figure  10  Stress  intensity  factors  at  the  tips  and  6i  of  two  unequal  cracks  in  an 
isotropic  FGM  under  fixed-grip  condition,  <jy y{x,0)  =  -E0e 0exp(ax),  k0  =  EqEq yjc[, 
Eo  =  E0/(l  -  i/2),  (c2  +  d)M  =  0.75,  i/  =  0.3. 
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Figure  1 1  Stress  intensity  factors  at  the  tips  ai  and  bi  of  two  unequal  cracks  in  an 
isotropic  FGM  under  fixed-grip  condition,  <ryy(x,  0)  =  -Eoeoexp(ax),  fa  =  EqEq^/c\, 
E0  =  Eo/(l  -  v2),  c2/ci  =  0.25,  v  =  0.3. 
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Figure  12  Stress  intensity  factors  at  the  tips  ax  and  bi  of  two  unequal  cracks  in  an 
isotropic  FGM  under  fixed-grip  condition,  ayy(x,  0)  =  -E0£oexp(ax),  fa  =  EoSoy/ci, 
Eq  =  Eq/{1  -  v2),  aci  =  0.5,  u  =  0.3. 
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ac\ 

Figure  13  Stress  intensity  factors  at  the  tips  aa  and  62  of  two  unequal  cracks  in  an 
isotropic  FGM  under  fixed-grip  condition,  cryy(x,0)  =  -Eoe0exp(ax),  fa  =  E0eoy/ci, 
Eq  =  Eq/(  1  -  v2),  02/ci  =  0.1,  v  =  0.3. 


(ci  +c2)/d 


Figure  14  Stress  intensity  factors  at  the  tips  a\  and  b\  of  two  unequal  cracks  in  an 
isotropic  FGM  under  fixed-load  condition,  ayy{x,  0)  =  -po,  fa  =  po^/ci",  oijc\  =  0.25, 
v  =  0.3. 
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Figure  15  Stress  intensity  factors  at  the  tips  a\  and  6iof  two  unequal  cracks  in  an  isotropic 
FGM under  fixed-load  condition,  ayy(x,0)  =  -po,  ko  -  poy/c[,  c2/ci  =  0.25,  u  =  0.3. 


OiC\ 

Figure  16  Stress  intensity  factors  at  the  tips  and  b2  of  two  unequal  cracks  in  an 
isotropic  FGM  under  fixed-load  condition,  cr^x,  0)  =  -po,  k0  =  Poyfci,  c2/c\  =0.1, 
v  =  0.3. 


ac\ 


Figure  17  Stress  intensity  factors  at  the  tips  ai  and  &iof  two  unequal  cracks  in  an 
orthotropic  FGM  under  fixed-grip  condition,  crw(x,0)  =  —  eoEo/6oexp(axi), 
ko  =  ( e0Eo/6o)y/c T,  c2/ci  =  0.25,  k  =  0.5,  v  =  0.3. 
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Figure  18  Stress  intensity  factors  at  the  tips  a\  and  £>iof  two  unequal  cracks  in  an 
orthotropic  FGM  under  fixed-grip  condition,  ayy(x,0)  =  -£QEo/6Qexp(axi), 
ko  =  {eoEo/6l)y/c{,  (c2  +  Ci)/d  =  0.75,  «  =  0.5,  v  =  0.3. 


aci 


Figure  19  Stress  intensity  factors  at  the  tips  ai  and  &iof  two  unequal  cracks  in  an 
orthotropic  FGM  under  fixed-grip  condition,  ayy(x,0)  =  -^o^o/5oexP(a:ci)> 
ko  =  (Eo£o/6o)\/ci>  (c2  +  ci )/d  =  0.5,  c2/ci  =  0.25,  v  —  0.3. 


(ci  +c2)/d 

Figure  20  Stress  intensity  factors  at  the  tips  aj  and  b\  of  two  unequal  cracks  in  an 
orthotropic  FGM  under  fixed-grip  condition,  cryy(x,0)  =  -Eo£o/^oexP(aa:i)> 
ko  =  (Eoeo/SDy/ci,  aci  =  0.5,  c2/ci  =  0.25,  u  =  0.3. 
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Figure  21  Stress  intensity  factors  at  the  tips  ax  and  bx  of  two  unequal  cracks  in  an 
orthotropic  FGM  under  fixed-grip  condition,  ayy(x,0)  = —E0£o/62exp(axi), 
ho  =  Eo£o/62y/ci,  ac\  =  0.5,  v  =  0.3,  k  =  0.5. 


x/a 


Figure  22  Crack  surface  displacement  for  an  isotropic  inhomogeneous  medium  under 
fixed-grip  condition,  disregarding  interpenetration  of  crack  surfaces, 
<ryy(x,0)  =  -E0ei(xi/a)exp(axi),  v0  =  2exa,  E0  =  E0/(  1  -  v2),  v  =  0.3. 
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Figure  23  Crack  surface  displacement  for  an  isotropic  inhomogeneous  medium  under 
fixed-grip  condition  in  the  presence  of  smooth  contact,  <ryy(x,  0)  =  -Eq£i  (x/a)exp(ax), 
vo  =  2eia,  E0  =  Eo/(l  -  v2),v  =  0.3,  £i  >  0. 


Figure  24  Crack  surface  displacement  for  an  isotropic  inhomogeneous  medium  under 
fixed-grip  condition  in  the  presence  of  smooth  contact,  ayy{x,  0)  =  -Eq€\  (x/a)exp(ax), 
Vo  =  2ei a,  Eq  =  Eq/(\  —  v2),v  =  0.3,  e\  <  0. 

35 


x/a 

Figure25  Crack  surface  displacements  for  an  infinite  isotropic  inhomogeneous  medium 
under  fixed-grip  condition  :  eyy  =  e0  +  £i(x/a),  e0M  =  P,  pc  =  0.4485,  aa  =  0.5, 
v  =  0.3. 


x/cx 

Figure  26  The  stress  distribution  in  the  plane  of  the  crack  in  an  inhomogeneous  isotropic 
medium  with  two  unequal  colinear  cracks  under  the  remote  loading  a™{x,  0)  =  e0EQeax, 
otc\  =  0.25,  c2/ci  =  0.1,  (ci  +  c2)/d  =  0.25,  cryy{x,  0):  the  total  stress,  a%y(x,  0)  :  the 
perturbation,  a^(x,  0):  stress  at  infinity  (the  applied  stress),  0yy  =  a%y  +  cr~. 
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Table  1  The  normalized  stress  intensity  factors  for  an  inhomogeneous  isotropic 
medium  under  crack  surface  tractions  given  by  Eq.  27  (  v  =  0.3  ). 


Plane  Strain 

fci(a) 

fci(-a) 

fci(a) 

fci(-a) 

ki  (a) 

fci(-a) 

OLCL 

Po\/a 

Poy/a 

Piy/a 

Piy/a 

P2y/a 

P2\/o 

P3y/a 

P3\/a 

0.00 

1.0 

1.0 

0.5 

-0.5 

0.5 

0.5 

0.375 

-0.375 

0.01 

1.0025 

0.9975 

0.5000 

-0.5000 

0.5006 

0.4994 

0.3750 

-0.3750 

0.10 

mwm'M 

0.9740 

0.4998 

-0.4998 

0.5060 

0.4935 

0.25 

1.0567 

0.9334 

0.4989 

-0.4986 

0.5142 

0.4833 

0.3745 

-0.3743 

0.50 

1.1062 

0.8667 

0.4962 

-0.4944 

0.5267 

0.4665 

0.3731 

-0.3722 

0.75 

1.1504 

0.8044 

0.4923 

-0.4878 

0.5379 

0.4505 

0.3711 

-0.3689 

1.00 

1.1902 

0.7480 

0.4876 

-0.4793 

0.5483 

0.4357 

0.3688 

1.50 

1.2598 

0.6536 

0.4771 

-0.4590 

0.5668 

0.4098 

0.3634 

-0.3542 

2.00 

1.3195 

0.5810 

0.4660 

-0.4371 

0.5833 

0.3881 

0.3580 

-0.3429 

Plane  Stress 

ki  (a) 

fci(-a) 

fci(a) 

ki(-a) 

k\(a) 

ki(-a) 

ki(a) 

ki(-a) 

Otd 

Po\fa 

Poy/a 

Piy/a 

Pi\fa 

P2\/a 

P2y/a 

Pzy/a 

Pzy/a 

0.00 

1.0 

l.o 

0.5 

-0.5 

0.5 

0.5 

0.375 

-0.375 

0.01 

1.0025 

0.9975 

0.5000 

-0.5000 

0.5006 

0.4994 

0.3750 

-0.3750 

0.10 

liEKM 

0.9737 

miLwim 

-0.4998 

0.5059 

0.4934 

0.25 

1.0553 

0.9324 

0.4989 

-0.4985 

-0.3720 

0.50 

1.1019 

0.8640 

0.4962 

0.3731 

-0.3720 

0.75 

1.1421 

0.8002 

-0.4870 

0.5359 

0.4494 

0.3712 

-0.3685 

Hi!il 

1.1774 

0.7427 

0.4343 

0.3689 

gik'telll 

1.2369 

0.6473 

0.4780 

0.3639 

|  2.00 

1.2862 

0.5746 

0.4680 

-0.4345 

0.5749 

0.3859 

0.3590 

-0.3415 

Table  2  The  normalized  stress  intensity  factors  for  an  inhomogeneous 
isotropic  medium  under  fixed  -  grip  condition  (  Eq  —  E0  for  plane  stress, 
Eq  =  Eq/(1-  v2)  for  plane  strain ,  v  =  0.3,  see  Eq.  30). 


i/  =  0.3  Plane  Strain  Plane  Stress 


aa 

fci(a) 

M-a) 

Ma) 

ki(-a) 

ki  (a) 

fci(-a) 

eoEoyfi 

toEoy/a 

CiEoy/a 

toEoy/a 

BEW m 

m mm 

0.00 

1.0 

1.0 

0.5 

1.0 

1.0 

0.5 

-0.5 

0.01 

1.0075 

0.9925 

0.5050 

1.0075 

0.9925 

0.5050 

-0.4950 

0.10 

1.0764 

0.9264 

0.5523 

-0.4523 

0.9261 

0.5523 

-0.4522 

0.25 

1.1986 

0.8230 

0.6402 

-0.3886 

0.8219 

0.6401 

-0.3885 

0.50 

1.4290 

0.6710 

0.8151 

-0.3010 

k mm 

0.6683 

0.8146 

-0.3009 

0.75 

1.7029 

0.5434 

-0.2325 

1.0311 

-0.2326 

1.00 

MikkVJi 

0.4379 

-0.1794 

1.2997 

-0.1795 

1.50 

0.2811 

2.0570 

-0.1064 

2.9012 

2.0482 

-0.1067 

2.00 

4.2933 

0.1786 

3.2297 

-0.0631 

4.2454 

0.1741 

3.2107 

-0.0633 
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Table  3  The  effect  of  Poison's  ratio  on  the  normalized  stress  intensity  factors 
for  an  inhomogeneous  isotropic  medium  (E0  —  Eq  for  plane  stress, 
Eo  =  Eo/(  1  -  v2)  for  plane  strain,  see  Eqs.  28  and  30  ). 


Fixed 


M-«) 


0.6615 


0.6641 


0.6672 


0.6710 


55 


0.6814 


Plane  Strain 


■Grip 


ki(a) 


k\(—a) 


ci  Ei 


-0.3009 


-0.3009 


Fixed-Load 


0.8131 


0.8137 


0.8143 


0.8151 


0.8162 


0.8176 


-0.3010 


-0.3011 


M«) 


€qEo\Jo. 


1.4129 


1.4170 


1.4208 


1.4245 


1.4280 


1.4314 


Fixed 


M-a) 


€qEq  y/a 


0.6615 


M-°) 


0.6724 


Plane  Stress 


•Grip 


M°) 


eiEo-s/a 


0.8131 


0.8136 


-0.3009 


-0.3009 


-0.3009 


-0.3010 


M«) 

fci(-a) 

ki(a) 

fcx(-a) 

Poy/a 

Piy/a 

Piy/a 

1.0906 

0.8570 

0.4961 

BigESSlI 

1.0950 

0.8597 

0.4961 

-0.4935 

1.1001 

0.4962 

Boembm 

1.1062 

0.8667 

0.4962 

-0.4944 

1.1137 

0.8715 

0.4963 

1.1234 

0.8775 

0.4964 

-0.4959 

Fixed-Load 

M®) 

fei(-a) 

M®) 

M-®) 

Poy/a 

Poy/a 

P\y/a 

P\yfa 

1.0906 

0.8570 

0.4961 

-0.4932 

1.0946 

0.8594 

0.4961 

B9BMa 

1.0983 

0.8618 

0.4961 

msEsi 

1.1019 

0.8640 

0.4962 

gjgrni 

1.1053 

0.8661 

0.4962 

-0.4943 

1.1085 

0.8682 

0.4962 

-0.4946 

Plane  Strain 

eta  =  1.0 

Fixed 

-Grip 

Fixed-Load 

V 

M«) 

*i  (-<*) 

M°) 

fci(-a) 

ki{a) 

fci(-a) 

M«) 

M~ ®) 

Poy/a 

Poy/a 

0.0 

1.9819 

0.4208 

1.2917 

-0.1798 

1.1438 

0.7291 

0.4888 

-0.4752 

0.1 

1.9963 

0.4256 

1.2948 

-0.1797 

1.1569 

0.7344 

0.4884 

-0.4763 

0.2 

2.0132 

0.4312 

1.2984 

-0.1795 

1.1722 

razTira 

0.4880 

-0.4777 

0.3 

2.0332 

0.4379 

1.3029 

-0.1794 

1.1902 

0.7480 

0.4876 

-0.4793 

0.4 

2.0576 

0.4459 

1.3084 

-0.1793 

0.4874 

-0.4813 

0.5 

2.0885 

0.4560 

1.3157 

-0.1793 

0.4873 

-0.4840 

Plane  Stress 

aa  —  1.0 

Fixed 

-Grip 

Fixed-Load 

V 

ki(a) 

fci(-a) 

fci(a) 

ki(-a) 

ki  (a) 

fci(-a) 

ki(a) 

ki(-a) 

P\y/a 

Ply/a 

0.0 

1.9819 

1.1438 

0.7291 

-0.4752 

0.1 

1.9949 

1.1556 

0.7339 

inrei 

0.2 

0.4292 

1.2971 

-0.1796 

0.3 

0.4331 

1.2997 

-0.1795 

0.4879 

-0.4781 

0.4 

2.0301 

0.4368 

1.3022 

-0.1794 

1.1874 

0.7468 

0.4877 

-0.4790 

0.5 

2.0407 

0.4404 

1.3046 

-0.1793 

0.4875 

-0.4799 
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Table  4  The  normalized  stress  intensity  factors  for  an  inhomogeneous  isotropic 
medium  under  fixed  -  load  condition  for  the  case  of  plane  strain  {v  =  0.3  ). 


Present 

Delale  and  Erdogan  [1983] 

aa 

ki{a) 

ki(-a) 

ki(a) 

ki(-a) 

ki(a) 

ki(-a) 

ki(a) 

fci(-a) 

Poy/a 

Po\/a 

Pi\fa 

Pi\fa 

P(*/a 

Po\/a 

Piy/a 

Pi\/ci 

0.00 

1.0 

1.0 

0.5 

-0.5 

1.0 

1.0 

0.5 

-0.5 

0.01 

1.0025 

0.9975 

0.5000 

agiKTiTtlBI 

1.003 

6.500 

mmm 

0.10 

1.0238 

0.9740 

0.4998 

1.026 

0.500 

-0.500 

0.25 

1.0567 

0.9334 

0.4989 

-0.4986 

1.061 

0.931 

0.498 

-0.499 

0.50 

1.1062 

0.8667 

0.4962 

-0.4944 

1.117 

0.863 

0.494 

-0.495 

0.75 

1.1504 

0.8044 

0.4923 

-0.4878 

1.170 

0.801 

0.489 

-0.489 

1.00 

1.1902 

0.4876 

-0.4793 

1.222 

0.745 

0.483 

Big  1:11 

Table  5  The  normalized  stress  intensity  factors  for  an  inhomogeneous  isotropic 
medium  under  fixed  -  load  condition  for  the  case  of  plane  stress  ( v  =  0.3  ). 


Present 

Delale  and  Erdogan  [1983] 

aa 

fci(a) 

ki(-a) 

ki(a) 

ki(-a) 

ki(a) 

53 

T— 1 

fci(a) 

ki(-a) 

PQy/a 

Poy/a 

Piy/a 

Pi\/cl 

P<y/a 

Poy/a 

Pi\/a 

Pi\/a 

0.00 

1.0 

1.0 

6.5 

1.0 

1.0 

0.5 

-0.5 

0.01 

1.0025 

0.9975 

0.5000 

0.997 

0.500 

-0.500 

0.10 

1.0235 

0.9737 

0.4998 

-0.4998 

0.973 

0.500 

-0.500 

0.25 

1.0553 

0.9324 

0.4989 

-0.4985 

1.060 

0.930 

0.498 

-0.499 

0.50 

1.1019 

0.8640 

-0.4941 

1.113 

0.861 

■E2EEI 

0.75 

1.1421 

0.8002 

0.4923 

-0.4870 

1.162 

0.797 

0.489 

-0.489 

1.00 

mwBiiim 

Wfliiwm 

-0.4781 

1.209 

0.740 

0.483 

Hfl! 

Table  6  The  normalized  stress  intensity  factors  for  an  inhomogeneous  isotropic 
medium  under  fixed  -  grip  condition  (  Eq  =  Eq  for  plane  stress, 
Eq  =  Eq/(1  —  i'2)  for  plane  strain  ). 


Plane  Strain 

Plane  Stress 

v  =  0.3 

Present 

Delale  and  Erdogan  [1983] 

Present 

Delale  and  Erdogan  [1983] 

aa 

*i(a) 

M-a) 

kt(a) 

M-«) 

ki(a) 

ki(-a)  | 

eoEoy/a 

eoEoy/a 

e0E0\/a 

eoEoy/a 

eoEo^/a 

eoEoy/a 

toEoy/a 

0.00 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

0.01 

1.0075 

0.9925 

1.008 

0.992 

1.0075 

0.9925 

1.008 

0.992 

0.10 

i 

KES9I 

EESSHI 

1.0761 

0.9261 

1.078 

0.925 

0.25 

POI 

KsEHiHi 

ntva 

HI 

KB9I 

0.50 

1.4290 

0.6710 

0.667 

0.75 

1.7029 

0.539 

0.5395 

1.713 

0.535 

1.00 

2.0332 

0.4379 

2.063 

2.048 

0.429 
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Table  7  Stress  intensity  factors  at  the  tips  x  =  a\  and  x  —  b\  of  two  unequal 
cracks  in  an  isotropic  FGM  for  the  case  of  plane  strain,  v  =  0.3, 
E0  =  E0/(l  -  it). 


Fixed-Grip 

<Tyy(x,  °)  =  -^oeoexp(ax),  ko  =  EoSo^/ci 

ac\  =  0.5 

oijc\  =  0.25 

02/ci  =  0.75 

o 

1—1 

II 

1—1 

o 

(ci  +  c2)/d 

MM 

ko 

MM 

ko 

MM 

ko 

MM 

ko 

MM 

*0 

MM 

ko 

0.0 

0.6710 

1.4290 

0.6710 

1.4290 

0.6710 

1.4290 

0.1 

0.6713 

1.4299 

0.6728 

1.4334 

0.6735 

1.4350 

0.2 

0.6720 

1.4319 

0.6766 

1.4435 

0.6791 

1.4495 

0.4 

0.6740 

1.4405 

0.6873 

1.4805 

0.6947 

1.5000 

0.6 

0.6773 

1.4634 

0.7030 

1.5601 

0.7171 

1.6032 

0.8 

0.6835 

1.5440 

0.7284 

1.7721 

0.7526 

1.8624 

0.975 

0.7034 

2.3181 

0.7883 

3.1421 

0.8314 

3.4198 

1.0 

0.7656 

oo 

0.9395 

OO 

1.0210 

OO 

Fixed-Load 

Gyy  PO,  k0  Po^/cT 

aci  =  0.5 

C2/C1  =  0.25 

c2/ci  =  0.75 

C2/C1  =  1.0 

(ci  +  C2)/d 

MM 

k0 

MM 

ko 

MM 

Aro 

MM 

ko 

MM 

ko 

MM 

kQ 

0.0 

0.8667 

1.1062 

0.8667 

1.1062 

0.8667 

1.1062 

0.1 

0.8667 

1.1062 

0.8667 

1.1062 

0.8667 

1.1062 

0.2 

0.8668 

1.1063 

0.8668 

1.1064 

0.8668 

1.1064 

0.4 

0.8674 

1.1088 

0.8687 

1.1124 

0.8688 

1.1125 

0.6 

0.8692 

1.1198 

0.8750 

1.1405 

0.8765 

1.1436 

0.8 

0.8735 

1.1686 

0.8893 

1.2445 

0.8945 

0.975 

0.8881 

1.7064 

0.9280 

2.0652 

0.9422 

2.1280 

1.0 

0.9334 

OO 

1.0258 

OO 

1.0503 

OO 

Table  8  The  normalized  stress  intensity  factors  for  two  unequal  cracks  in  an 
infinite  inhomogeneous  orthotropic  medium  under  uniform  crack  surface 
pressure,0’yy(x,O)  =  —  po,E(xi)  =  Eoex-p(cxxi),u  =  0.3,  (ci  +  c2 )/d  =  0.5. 


k  =  —0.25,  C2/C1  ~  0.25 

k  =  0.0,  cs/ci  =  0.25 

a  ci 

MM 

P0y/ci 

MM 

POy/cl 

MM 

PQy/ci 

MM 

POyfcl 

aci 

MM 

POy/Cl 

MM 

Po^/cl 

MM 

POyfcl 

MM 

POy/ci 

0.0 

1.00429 

1.01010 

0.55216 

0.54146 

0.0 

1.00429 

1.01010 

0.55216 

0.54146 

0.01 

1.00171 

1.01253 

0.55240 

0.54232 

0.01 

1.00171 

1.01252 

0.55240 

0.54231 

0.10 

0.55437 

0.54991 

0.10 

liECTESI 

1.03357 

0.55421 

0.54976 

0.25 

0.55710 

0.56227 

0.55646 

0.56166“ 

0.50 

0.87545 

1.12420 

0.56016 

0.87163 

1.11809 

0.55856 

0.58038 

0.75 

0.81612 

1.17694 

0.56131 

0.60011 

0.75 

0.81038 

1.16556 

0.55884 

0.59774 

LOO 

0.76203 

1.22710 

0.56061 

0.61679 

1.00 

0.75489 

1.20973 

0.55749 

0.61377 

1.25 

0.55823 

0.63194 

1.25 

0.70538 

1.25104 

0.55471 

0.62849 

1.50 

0.55437 

0.64562 

rail 

0.66160 

1.28987 

0.55069 

0.64197 

1.75 

1.36367 

0.54923 

0.65794 

ita 

0.62305 

1.32655 

0.54562 

0.65430 

2.00 

0.59792 

1.40519 

0.54304 

0.66901 

1.36137 

0.53969 

0.66557 

Table  8  (Continued  ) 


K  =  0.5,  C2  /c\ 

=  0.25 

K  = 

L.O,  oil  ci  =  0.25 

Mai) 

MM 

Ma2) 

fcl(&2) 

Mai) 

MM 

Ma2) 

MM 

CtC\ 

PO-y/cT 

poy/cT 

POy/cT 

Po\/ci 

OtC\ 

P oy/ci 

PO-y/ci 

Poi/ci 

Paifcl 

0.0 

1.00429 

1.01010 

0.55216 

0.54146 

0.0 

1.01010 

0.55216 

0.54146 

0.01 

1.00170 

1.01251 

0.55239 

0.54231 

0.01 

1.00170 

1.01251 

0.55239 

0.54231 

0.97746 

1.03312 

0.55404 

0.54960 

0.10 

0.97725 

1.03289 

0.55395 

0  .54952 

0.25 

0.93578 

1.06459 

0.55575 

0.56099 

0.25 

0.93491 

1.06348 

0.55536 

0.56061 

0.50 

0.86754 

1.11153 

0.55669 

0.57857 

0.50 

0.86538 

1.10805 

0.55563 

0.57753 

0.75 

0.80420 

1.15325 

0.55579 

0.59474 

0.75 

0.80091 

1.14668 

0.55399 

0.59294 

1.00 

0.74718 

1.19085 

0.55341 

0.60968 

1.00 

0.74307 

1.18072 

0.55090 

1.25 

IflRSTEEl 

1.22509 

0.54981 

0.62349 

1.25 

0.69204 

1.21113 

0.54669 

0.62021 

1.50 

0.65236 

1.25657 

0.54521 

0.63626 

1.50 

0.64742 

1.23862 

0.54157 

0.63234 

1.75 

0.61357 

1.28579 

0.53981 

0.64809 

1.75 

0.60850 

1.26377 

iinaa 

0.64362 

2.00 

0.57964 

1.31314 

0.53375 

0.65906 

MfltI 

0.57457 

1.28702 

0.52942 

0.65412 

K  = 

5.0,  c2/< 

:i  =  0.25 

_ 

OtC\ 

Mai) 

MM 

ki  (i a 2) 

MM 

POy/cT 

P0y/ci 

P0-y/c7 

POy/cT 

0.0 

1.00429 

1.01010 

0.55216 

0.54146 

0.01 

1.00169 

1.01250 

0.55239 

0.54231 

0.10 

0.97681 

1.03240 

0.55376 

0.54933 

0.25 

0.93307 

1.06111 

0.55451 

0.55979 

0.50 

0.86075 

1.10060 

0.55317 

0.57510 

0.75 

0.79385 

1.13250 

0.54967 

0.58855 

1.00 

0.73421 

0.60059 

1.25 

0.68203 

1.18078 

0.53860 

0.61150 

1.50 

0.63676 

1.19948 

0.53179 

0.62148 

1.75 

0.59757 

1.21562 

0.52447 

0.63068 

1.22978 

0.51680 

0.63921 

C2/C1  =  1.0 


MM 

MM 

ki  (02) 

PO-y/cT 

PO-y/cT 

PO-y/cT 

1.02796 

1.04795 

1.04795 

1.02477 

1.04978 

1.0461 

.99602 

1.06556 

1.0288 

.94961 

1.09128 

0.99839 

.87927 

1.13545 

0.94573 

.81729 

1.18161 

0.89313 

.76231 

1.22847 

0.84244 

.71346 

1.27474 

0.79476 

.67013 

1.31969 

0.75064 

.63179 

1.36299 

0.71025 

.59791 

1.40459 

0.67350 

MM 


POl/oi 


1.02796 


1.0311 


1.0591 


1.1045 


1.1764 


1.2432 


1.3046 


1.3608 


1.4121 


1.5024 


0.0,  C2/C1  =  1.0 


Mai)  MM 


Poy/ci  PoJcl 


1.02796  1.04795 


1.02476  1.04977 


IHCSIEEMIiTSEEll 


0.94786 


0.87537 


0.75  0.81153 


1.00  0.75517 


1.25  0.70543 


1.50  0.66159 


0.58914 


1.08894 


1.12905 


1.17010 


1.21108 


1.25099 


1.28932 


1.32590 


1.36080 


1.0276  I  1.0579 


k  —  0.5,  c2/ci  =  1.0 


Ma  2)  M&  2)  Mgi)  fci  (bx )  Ma2)  MM 

Po-y/cT  Poy/ci  Po-y/cT  Poy/ci  Po-y/cT  Po-y/al 


1.04795  1.02796  0.0  1.02796  1.04795  1.04795  1.02796 


1.0461  1.0311  0.01  1.02475  1.04975 


1.0282  1.0585  0.10  0.99499  1.06435 


0.99606  1.0190?  0.25  0.94598  1.08641 


0.94040  1.1691  0.50  0.87115  1.12210 


0.88560  1.2306  0.75  0.80530  1.15758 


0.83370  1.2865  1.00  0.74746  1.19212 


0.78567  1.3373  1.25  0.69673  1.22504  0.77421  1.3092 


0.74185  1.3835  1.50  0.65235  1.25606  0.73029  1.3489 


0.70217  1.4255  1.75  0.61356  1.28518  0.69109  1.3848 


0.66641  1.4640  2.00  0.57964 


0.93434  1.1608 
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K  =  1.0,  02 

Id  =  1.0 

K  S5S  5.0,  02/ci  =  1.0 

aci 

MM 

POy/ci 

*1(6!) 

flJ-y/cT 

Ma  2) 
Po-Jci 

*1(62) 

Poy/ci 

OLC\ 

MM 

Poy^ 

*1(61) 

Poy/ci 

*1(02) 

Poy^ 

*1(62) 

Poy/ci 

0.0 

1.02796 

1.04795 

1.04795 

1.02796 

0.0 

1.02796 

1.04795 

1.04795 

1 .02796 

0.01 

1.02474 

1.04975 

1.0461 

■RricW 

0.01 

1.02473 

1.04973 

1.0461 

1.0311 

0.10 

0.99471 

1.06402 

1.0272 

1.0575 

0.10 

0.99412 

1.06333 

1.0265 

1.0568 

0.25 

0.94499 

1.08506 

0.99213 

1.0974 

0.25 

0.94284 

1.08218 

0.98913 

1.0940 

0.50 

0.86890 

1.11838 

0.93097 

1.1562 

0.50 

liTilgTCT 

1.11032 

0.92338 

1.1457 

0.75 

0.80198 

1.15086 

0.87168 

1.2075 

0.75 

0.79481 

1.13628 

0.85991 

1.1881 

1.00 

0.74334 

1.18194 

0.81682 

1.2527 

1.00 

0.73445 

1.15979 

0.80183 

1.2234 

1.25 

0.69209 

1.21107 

0.76727 

1.2927 

1.25 

0.68207 

1.18065 

0.75005 

usual 

1.50 

0.64741 

1.23812 

0.72308 

1.3285 

1.50 

0.63675 

1.19898 

0.70449 

1.2784 

1.75 

0.60849 

1.26318 

0.68390 

1.3606 

1.75 

0.59756 

1.21507 

0.66462 

1.3002 

2.00 

0.57456 

1.28649 

0.64923 

1.3897 

2.00 

0.56360 

1.22929 

0.62976 

1.3194 

Table  9  The  normalized  stress  intensity  factors  for  two  unequal  cracks  in  an 
infinite  inhomogeneous  orthotropic  medium  under  uniform  crack  surface 
pressure,  ayy(x,0)  =  —  po,  E(x i)  =  EoQxp(axi),  v  =  0.3,  k  =  0.5. 


02/ci  =  0.25,  (ci  -f  c^/d  =  0.75 

aci 

*1(01) 

P0y/Cl 

MM 

Poi/c7 

ki(a2) 

Poy/ci 

MM 

POy/ci. 

0.0 

0.66867 

0.60894 

0.01 

1.00880 

1.05090 

0.66935 

0.61018 

1.07123 

0.67479 

0.62088 

0.25 

0.94123 

1.10193 

0.68212 

0.63757 

0.50 

0.87140 

1.14689 

0.69093 

0.66332 

0.75 

0.80678 

0.69652 

0.68704 

1.00 

0.74884 

1.22084 

0.69956 

0.70907 

1.25 

0.69772 

■resnai 

0.70055 

0.72964 

1.50 

0.65298 

1.28079 

0.69990 

liirraaf 

1.75 

0.61394 

1.30729 

0.69792 

0.76707 

2.00 

0.57986 

1.33211 

0.78420 

C2/C1  =  0.75,  (ci  4*  Oi)/d  =  0.75 

aci 

*1  (ai ) 
P0y/ci 

MM 

Pov/ci 

MM 

POy/ci 

*1(62) 

PD-v/ci 

0.0 

1.04964 

1.14076 

1.04573 

0.94551 

0.01 

1.04636 

1.14252 

1.0452 

0.94823 

0.10 

1.01574 

1.15652 

1.0389 

0.97155 

0.25 

0.96426 

1.17566 

1.0244 

1.0077 

0.50 

0.88389 

1.20191 

0.99440 

1.0632 

0.75 

0.81314 

1.22513 

0.96126 

1.1138 

1.00 

1.24726 

0.92738 

1.1604 

1.25 

0.89413 

1.2034 

U*il 

0.65361 

1.29073 

0.86230 

1.2432 

1.75 

0.61421 

WSFETil 

0.83227 

1.2802 

2.00 

0.57996 

1.33373 

0.80422 

1.3146 

C2/C1 

=  1.0, 

(ci  +  c2)/d  =  0.75 

aci 

MM 

*1(61) 

*1  (M 

MM 

pov/cj 

P0y/Cl 

P0y/Cl 

P0y/Cl 

0.0 

1.06844 

1.17318 

1.17318 

1.06844 

0.01 

1.06475 

1.17451 

1.1717 

1.0720 

0.10 

1.03045 

1.18417 

1.1565 

1.1028 

0.25 

0.97370 

1.19589 

1.1262 

1.1500 

0.50 

0.88781 

1.21176 

1.0707 

1.2213 

0.75 

0.81450 

1.22790 

1.0150 

1.2853 

1.00 

0.75226 

1.24581 

0.96240 

1.3434 

1.25 

0.69915 

1.26541 

0.91427 

1.3963 

1.50 

0.65354 

1.28619 

0.87088 

1.4448 

1.75 

0.61414 

1.30762 

0.83206 

1.4894 

2.00 

0.57991 

1.32931 

0.79742 

1.5307 
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02/ci  =  0.25,  aci  =  0.5 

C2/C1  =  0.75 ,  aci  =  0.5 

C1+C2 

d 

fci(ai) 

POy/cl 

MM 

POy/ci 

MM 

Poy/ci 

MM 

Poy/ci 

Cl  +  C2 

d 

MM 

POy/ci 

MM 

Po-s/ci 

MM 

POy/ci 

MM 

P0y/ci 

0.0 

0.8661 

1.1053 

0.4331 

0.5527 

0.0 

0.8661 

1.1053 

0.7501 

0.9572 

0.1 

0.86614 

1.10527 

0.48805 

0.51930 

0.1 

0.86614 

1.10526 

0.78301 

0.94138 

0.2 

0.86619 

1.10539 

0.49731 

0.52839 

0.2 

0.86622 

1.10545 

0.79236 

0.95090 

0.3 

0.86640 

1.10607 

0.51092 

0.54096 

0.3 

0.86674 

1.10694 

0.80597 

0.96313 

0.4 

0.86683 

1.10783 

0.53004 

0.55743 

0.4 

0.86809 

1.11143 

0.82517 

0.97818 

0.5 

0.86754 

1.11153 

0.55669 

0.57857 

0.5 

0.87055 

1.12107 

0.85224 

0.99643 

0.6 

0.86863 

1.11885 

0.59451 

0.60569 

0.6 

0.87439 

1.13943 

0.89137 

1.0187 

0.7 

0.87026 

1.13374 

0.65093 

0.64115 

0.7 

■ana 

1.17380 

0.95114 

1.0464 

0.8 

0.87287 

1.16752 

0.74431 

0.68972 

0.8 

0.88862 

1.24313 

1.0531 

1.0828 

0.9 

0.87776 

1.26963 

0.94244 

0.76440 

0.9 

0.90295 

1.42220 

1.2792 

1.1376 

0.95 

0.88273 

1.44262 

1.1969 

0.82657 

0.95 

0.91588 

1.68871 

1.5825 

1.1836 

1.0 

0.9326 

OO 

OO 

1.2606 

1.0 

1.0244 

OO 

OO 

1.5458 

02/01  =  1.0,  aci  =  0.5 

M«i) 

PD-v/q 

MM 

MM 

po^c\ 

MM 

Po^c. I 

0.8661 

1.1053 

0.8661 

1.1053 

0.86614 

0.86986 

1.1096 

0.86620 

1.10541 

0.87854 

1.1187 

0.86671 

1.10682 

0.89116 

1.1302 

0.86823 

1.11153 

0.90900 

1.1441 

0.93434 

1.1608 

0.87588 

0.97136 

1.1810 

IBEEiKil 

1.18072 

1.0287 

1.2061 

0.89377 

1.25695 

1.1282 

1.2392 

0.91169 

1.45009 

1.3533 

1.2894 

1.6602 

1.3323 

1.0561 

OO 

OO 

1.6792 

Table  10  Stress  intensity  factors  and  the  contact  zone  sizes  for  an  isotropic 
inhomogeneous  medium  under  fixed  grip  remote  bending,  v  =  0.3, 
kQ  —  E0£iy/a ,  cryy(x,0)  =  -E0£i(x/a)exp(ax),E0  =  E0/(l  ~  v2)- 


£l  >  0 

£1  <  0 

No  contact 

ota 

6/a 

M°) 

6/a 

fcj(-a) 

Ma) 

ki(-a) 

—Eo£i  -y/a 

Eoe,^ 

E0eis/a 

0.0 

0.3333 

0.5443 

0.3333 

0.5443 

0.5 

-0.5 

iiy 

0.3405 

0.5964 

0.3258 

0.4965 

0.5523 

-0.4523 

0.25 

0.3508 

0.6832 

liKlEIil 

0.4320 

0.6402 

-0.3886 

0.5 

wiKmm 

0.8554 

0.3419 

0.8151 

-0.3010 

1.0 

1.3359 

0.2585 

0.2130 

1.3029 

-0.1794 

WEM 

0.4234 

2.0827 

0.2269 

0.1320 

2.0570 

-0.1064 

2.0 

0.4476 

3.2488 

0.2000 

0.0815 

3.2297 

-0.0631 
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Table  1 1  Stress  intensity  factors  and  contact  zone  sizes  for  an  isotropic  inhomogeneous 
medium  under  remote  strain  eyy  =  eo  +  £1  {x/a)  =  £\{P  +  x/a), (Fig.  25),  a  =  0.5, 
V  =  0.3,  £q/£\  =  /?,  £\  >  0. 


Si 

b/a 

&i(a) 

k\  (-a) 

Efciy/i 

E$e\y/a 

0.0 

0.3669 

0.8554 

0.0 

0.25 

0.7184 

1.1798 

0.0 

0.4485 

1 

1.4561 

0.0 

0.45 

1 

1.4582 

0.0010 

0.50 

1 

1.5297 

0.0345 

0.625 

1 

1.7083 

0.1184 

0.75 

1 

1.8869 

0.2023 

1.00 

1 

2.2442 

0.3700 
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